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Objectives 

This  project  was  set  to,  (a)  investigate  the  Kinematics  of  damage  before  crack 
initiation,  (b)  formulate  kinetic  equations  for  damage  growth,  and  (c)  derive 
criteria  for  crack  initiation  in  an  amorphous  material.  From  the  research  so 
far,  the  following  Important  results  have  emerged: 

1 .  Damage  accompanies  the  process  of  crack  initiation.  The  damage  zone 
consist  of  a  core  zone  with  highly  localized  damage  and  a  peripheral  less 
dense  damage  zone. 

2.  While  the  pattern  of  the  peripheral  zone  is  depended  upon  the  loading 
history,  the  density  in  the  core  zone  is  independent  of  the  loading  history. 

3.  The  kinematics  of  damage  before  crack  initiation  can  be  described  by  an 
affine  transformation  of  the  space  variables.  Accordingly,  damage  evolution 
can  be  described  by  the  translation,  isotropic  expansion  and  homogeneous 
distortion  of  the  damage  zone. 

4.  The  time  to  crack  initiation  is  related  to  the  stress  level  with  an 
exponentially  decaying  relationship.  The  results  agree  with  fracture  models 
based  on  kinetic  theories  for  submicoscopic  flaw  nucleation. 

The  results  of  the  present  studies  and  those  reported  in  [17,18]  indicate 
certain  similarities  between  the  characteristics  of  damage  before  crack  initiation 
and  during  slow  crack  growth.  These  are,  (a)  damage  evolution  can  be 
described  by  a  linear  transformation,  (b)  the  densities  of  damage  within  a  core 
zone  at  crack  initiation  and  during  slow  crack  growth  are  of  the  same  level. 

Relevant  publications 

The  work  so  far  has  resulted  in  two  papers  that  have  been  submitted  for 
publication: 

1 .  J.  Botsis  and  C.  Huang,  'Experimental  and  Analytical  Studies  in  Crack 
Initiation  under  Fatigue:  Part  I  -  Experimental',  International  Journal  of  Fracture, 
Submitted,  August,  1993. 

2.  J.  Botsis  and  C.  Huang,  'Experimental  and  Analytical  Studies  in  Crack 
Initiation  under  Fatigue:  Part  II  -  Analysis',  International  Journal  of  Fracture, 
Submitted,  August,  1993. 

Future  Work 

A.  Energy  calculations  due  to  damage 

B.  Kinetic  equations  for  damage  growth 

C .  Criteria  fro  crack  initiation 


The  results  obtained  so  far  are  given  in  the  subsequent  pages. 
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Summary 


Mechanistic  investigations  of  damage  evolution  before  crack  initiation  in  an 
amorphous  polymer  show  that  damage  consists  of  a  core  of  highly  dense 
crazing  and  a  peripheral  less  dense  zone  of  crazing.  Damage  characterization 
is  carried  out  at  consecutive  configurations  of  the  damage  zone.  Analysis  of  the 
kinematics  of  damage  at  different  times  involves  comparisons  of  the  inertia 
moments  of  damage  distributions.  The  results  indicate  that  damage  evolution 
between  consecutive  configurations  can  be  approximated  by  a  linear 
transformation  of  the  space  variables.  Thus,  the  process  of  damage  growth  can 
be  described  by  translation  and  deformation  of  the  damage  zone.  The  growth 
rates  of  the  damage  zone  movements  decrease  until  crack  initiation.  In  all 
cases,  the  average  densities  exhibit  a  damping  type  behavior  with  the  number 
of  cycles.  The  crack  initiates  within  a  core  zone  immediately  ahead  of  the  stress 
concentrator.  The  experimental  results  suggest  that  damage  density  within  the 
core  zone  is  independent  of  the  loading  conditions  considered  herein.  This 
value  is  approximately  equal  to  the  damage  density  around  the  crack  tip  during 
slow  crack  propagation.  The  crack  length  at  initiation  is  found  to  increase 
exponentially  with  the  stress  level.  A  simple  decaying  exponential  relationship 
relates  the  crack  initiation  times  and  the  applied  stress  level.  This  result  is 
consistent  with  the  fracture  models  based  on  absolute  reaction  theories. 
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1.  Introduction  and  Background 

Prediction  of  life  time  of  a  structural  component  subjected  to  fatigue  or  creep 
loads  assumes  a  position  of  prominer  ^e  in  engineering  design.  It  becomes 
even  more  important  when  considering  the  continuous  increase  in  complexity  of 
engineering  structures  and  the  large  costs  resulting  from  premature  failures. 
Thus,  characterization  of  the  mechanisms  leading  to  crack  initiation  and  growth 
is  an  issue  of  central  importance  to  engineers  and  material  scientists. 

Contemporary  approaches  to  life  time  prediction  consider  a  fracture 
process  in  three  phases:  crack  initiation,  slow  crack  growth,  and  fast  fracture. 
The  latter  phase  is  very  short  in  duration  and  thus,  the  time  of  crack  initiation 
and  that  of  slow  crack  growth  account  for  the  useful  time  of  a  structure. 

Significant  efforts  have  been  devoted  to  studying  slow  crack  propagation  in 
various  materials.  This  is  demonstrated  by  the  vast  amount  of  experimental  data 

and  several  models  that  have  been  proposed  to  describe  crack  growth  with  the 
stress  intensity  factor  Kv  or  energy  release  rate  J1  [1 ,  2j.  These  analyses  play 

an  important  role  in  life  time  prediction.  However,  experimental  and  theoretical 
studies  of  crack  initiation  have  received  limited  attention  in  the  literature. 
Moreover,  fracture  mechanics  approaches  to  fatigue  do  not  consider  crack 

initiation  [1].  Instead,  the  concept  of  crack  growth  threshold  is  employed  [3,  4]. 
Values  of  fatigue  threshold  indicate  the  cyclic  stress  intensity  factor  AKth,  below 

which  long  cracks  remain  dormant.  The  corresponding  AKth  is  taken  as  a 

material  parameter.  Although  the  study  of  fatigue  thresholds  and  near  threshold 
behavior  has  highlighted  the  effects  of  load  history  in  influencing  crack  growth, 
and  has  provided  insight  into  the  different  mechanisms  of  crack  closure,  the 
threshold  idea  has  met  with  limited  success  in  engineering  design.  This  stems 
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from  the  fact  that  the  use  of  a  fatigue  threshold  represents  a  conservative  design 
criterion  and  from  questions  regarding  as  a  material  parameter  [5-7]. 

Mechanistic  investigations  in  various  materials  have  shown  that  fatigue 
crack  initiation  and  growth  are  preceded  by  various  forms  of  damage  that 
nucleate  at  heterogeneities  within  the  material.  In  metals  initiation  is  related  to 
intense  slip  processes,  extrusions-intrusions  and  persistent  slip  bands  [8,  9]. 
Microcracks  nucleated  at  inclusions  or  second  phase  particles  in  the  case  of 
smooth  specimens  have  also  been  reported  by  several  researchers. 
Observations  of  the  same  process  in  both  notched  and  smooth  specimens  of 
polymers  have  emphasized  the  role  of  crazing  in  amorphous  and  semi  • 
crystalline  polymers  [10-13].  Damage  nucleation  and  growth,  whether  in  the 
form  of  microcracks,  voids,  homogeneous  transformations,  crazes,  etc.,  are 
processes  that  absorb  energy,  energy  that  otherwise  would  be  available  to 
drive  the  crack.  Therefore  the  nature  and  extent  of  damage  determine  the  time 
to  crack  growth  and  toughness  of  the  material.  Under  these  circumstances,  use 

of  fracture  mechanics  parameters  may  not  be  always  justified.  That  is,  when  the 
restrictions  of  ^-dominance'  or  ,J1 -dominance*  are  not  met,  K1  and  J1  fail  to 

correlate  with  the  fracture  process.  Examples  are  large  scale  plasticity, 
microcrack  interaction,  the  behavior  of  small  cracks,  etc. 

Recognizing  the  importance  of  damage  on  crack  initiation  and  growth, 
several  analytical  approaches  have  been  developed  to  account  for  its  effects  on 
stress  intensity  factor  and  energy  release  rates  (for  a  recent  review  see  [1 4]  and 
references  therein).  It  is  difficult,  however,  to  assess  the  contributions  of  crack 
damage  interactions,  describe  the  local  kinetics,  and  the  effects  of  the 
microstructure  when  a  large  number  of  microdefects  is  present.  These  powerful 
mechanisms  are  the  main  obstacles  in  modeling  fracture  in  a  large  class  of 
modern  engineering  materials. 


One  way  to  treat  fracture  phenomena  where  dense  damage  accompanies 
crack  initiation  and  growth  is  to  use  the  ideas  of  self  -  similarity.  Drawing  on  the 
success  of  statistical  self  -  similarity  of  vortex  cascade  in  understanding 
turbulence,  it  is  proposed  that  the  same  principles  can  be  applied  to  damage 
and  fracture  of  engineering  materials  [15,16].  Indeed,  experimental 
observations  in  a  number  of  materials  have  shown  that  damage  growth  displays 
certain  invariant  characteristics.  For  example,  the  evolution  of  a  normalized 
pore  size  distribution  in  polycrystalline  steel  under  creep  follows  a  pattern  that  is 
independent  of  time  and  strain  level  [15].  Furthermore,  it  has  been  shown  that 
damage  growth  within  a  process  zone  during  fatigue  fracture  can  be  described 
by  a  self  -  similar  transformation  of  the  space  variables  [17].  Modeling  of  such 
fracture  processes  has  been  made  by  introducing  additional  kinematic 
parameters  on  the  basis  of  self  *  similarity  [16]. 

Formulations  of  constitutive  equations  for  damage  growth  and  criteria  for 
crack  initiation  are  well  recognized  by  now  to  be  important  problems  in  fracture 
research.  Although,  important  progress  has  been  made  in  this  direction,  a 
number  of  questions  still  remain  to  be  properly  addressed  and  resolved  for  a 
better  understanding  of  fracture.  In  this  project,  inquiries  into  the  characteristics 
of  damage  evolution  before  crack  initiation  and  criteria  for  crack  initiation  are 
addressed.  In  the  first  phase,  experimental  studies  on  damage  evolution  are 
reported  using  amorphous  polystyrene  as  a  model  material.  To  facilitate 
experimental  observations,  the  site  of  crack  initiation  is  located  in  space  by 
inducing  a  60°  -V  notch  onto  the  mid-span  of  the  specimen  edge.  Hence, 
efforts  are  concentrated  on  characterizing  damage  dissemination,  the 
associated  kinematic  parameters,  and  the  time  to  crack  initiation.  In  the  second 
phase,  the  energy  release  rate  due  to  damage  growth  is  evaluated  using  of  a 
semi  -  empirical  method  and  experimental  measurements.  On  the  basis  of  the 


9 


experimental  results  damage  growth  in  front  of  the  notch  tip  is  modeled  with  a 
first  order  reaction  equation. 

2.  Experimental  Procedures 
2.1  Material  and  Specimens 

Commercially  available  amorphous  polystyrene  (PS)  from  Transilwrap, 
Chicago,  IL,  was  employed  in  the  studies  of  crack  initiation.  The  material  was 
received  in  the  form  of  sheets  with  dimensions  of  200x250  mm.  Strips  of 
150x22x0.20  mm  were  cut  from  these  sheets  with  a  razor  blade  and 
sandwiched  between  two  pieces  of  acrylic  with  dimensions  150x21x6.5  mm. 
The  specimens  were  squared  up  to  1 50x20x0.20  mm  with  the  use  of  a  fly  cutter 
on  a  milling  machine.  The  speeds  of  the  fly  cutter  and  the  longitudinal  power 
feed  of  the  working  table  that  supports  the  specimens  were  11  rotation-sec1, 
and  0.23  mm-sec*1,  respectively.  First,  sections  of  about  0.25  mm  thick  were  cut 
off  the  block  of  specimens  until  the  edges  of  the  specimens  were  on  the  same 
plane.  The  final  few  cuts  were  in  the  order  of  0.01  mm  so  that  the  induced 
damage  was  minimized.  The  block  was  then  flipped  over  and  the  same 
procedure  was  repeated  until  the  desired  width  was  achieved.  The  machined 
edges  were  metallographically  polished  to  a  0.5pm  finish  to  prevent  formation 
of  edge  crazes.  A  double  angle  cutter  of  60°  was  used  to  notch  the  specimens 
prior  to  their  removal  from  the  block.  Subsequently  the  specimens  were 
removed  from  the  block  and  washed  carefully  with  distilled  water.  Note  that  this 
procedure  of  specimen  preparation  ensures  identical  notch  tip  geometries. 
Finally,  the  specimens  were  annealed  in  a  temperature  of  10°C  degrees  lower 
than  Tg  (glass  transition  temperature)  for  48  hours  and  then  allowed  to  slowly 

cool  down  to  room  temperature.  Annealing  was  aimed  at  relieving  any  residual 
stresses  and  healing  any  damage  formed  during  the  cutting  procedures. 
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2.2  Experimental  Methods 

Tension  -  tension  fatigue  experiments  were  conducted  on  an  instron  Testing 
System  in  laboratory  environment  at  ambient  temperature.  All  experiments 
were  performed  under  load  controlled  mode  with  sinusoidal  wave  form.  The 
evolution  of  damage  around  the  notch  tip  was  observed  by  means  of  a  traveling 
optical  microscope  attached  to  the  Instron  Testing  System.  The  fracture 
process  was  recorded  using  a  motor  driven  camera  that  was  attached  to  the 
microscope.  Craze  distribution  was  evaluated  from  optical  micrographs  of 
sectioned  specimens  (approximately  10*20^m  thick)  which  were  prepared  by 
standard  metallographic  and  polishing  procedures.  Damage  density 
distributions  were  obtained  by  covering  the  micrographs  of  the  polished 
sections  with  a  mesh  of  rectangles.  The  size  of  a  typical  mesh  depended  upon 
the  extent  of  damage  around  the  notch  tip  and  was  about  N1xN2»  10x15 
(Figure  1 ).  In  all  cases  the  dimensions  of  the  rectangles  were  approximately 
12x1 2pm.  In  each  rectangle,  the  number  of  crazes  was  counted.  Craze  density 

was  evaluated  as  number  per  unit  area  p  [#/mm2]  or  as  p  [mm2/mm3].  p 

abt 

represents  the  amount  of  area  of  craze  mid  planes  per  unit  volume,  n  is  the 
number  of  intersections  of  crazes  with  the  vertical  test  line  at  the  respective 
rectangle,  a  and  b  were  the  height  and  width  of  a  rectangle,  respectively,  and  t 
is  the  specimen  thickness. 

The  most  direct  way  to  experimentally  measure  damage  growth  rates  are 
by  examination  of  sectioned  samples  that  have  been  exposed  to  a  well 
controlled  loading  history.  In  this  way  damage  growth  is  "frozen  in"  at  different 
degrees  of  development.  Thus,  to  obtain  damage  distribution  and  growth  rates 
before  initiation,  the  experiments  designed  for  damage  evolution  studies  were 
interrupted  at  appropriate  time  intervals.  Subsequently,  the  specimens  were 
polished  and  the  number  of  crazes  in  each  rectangle  was  evaluated. 
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The  approach  of  characterizing  damage  described  in  the  preceding 
paragraphs  implies  that,  at  a  given  time,  damage  distribution  around  the  stress 
concentrator  is  reproducible  for  a  certain  set  of  loading  conditions  and 
specimen  geometry.  Damage  accumulation,  however,  is  a  stochastic  process. 
Accordingly,  damage  distribution  should  result  from  averages  of  a  sufficiently 
large  number  of  distributions  observed  on  identical  specimens.  However,  such 
an  effort  is  experimentally  difficult  to  carry  out.  Commercially  available  PS  is  a 
relatively  'ductile'  material.  This  results  in  a  large  extent  of  crazing  during  crack 
initiation  and  slow  crack  growth  [17,18].  Thus,  damage  within  a  zone  in  this 
material  is  amenable  to  a  statistical  treatment.  To  examine  the  reproducibility  of 
damage  distribution,  two  experiments  were  performed  under  the  same  loading 
conditions  and  specimen  geometry. 

The  time  to  crack  initiation  was  examined  under  the  loading  histories 
shown  in  Table  I.  Crack  initiation  was  considered  when  a  crack,  however  small, 
would  appear  at  the  notch  tip.  This  event  was  observed  with  the  use  of  a  Nikon 
optical  microscope  under  a  magnification  of  x200.  It  has  been  shown  that  crack 
initiation  time  in  PS  under  fatigue  loads  exhibits  a  variability  of  about  10%  [19]. 
This  scatter  was  considered  small  and  thus,  efforts  were  concentrated  at 
investigating  damage  growth  before  initiation  under  three  loading  conditions, 
which  are  referred  to  as  Cl ,  C2,  and  C3.  Four  experiments  were  conducted 
under  Cl  and  three  under  each  of  C2  and  C3  conditions.  The  loads  and 
frequency  for  these  three  sets  of  experiments  are  shown  in  Table  I. 

Due  to  difficulties  associated  with  measuring  an  accurate  crack  size  at 
initiation  using  an  optical  technique,  the  tests  were  interrupted  shortly  after 
crack  initiation.  Subsequently,  the  specimens  were  pulled  to  fracture.  A  crack 
size  at  initiation  was  assessed  by  examination  of  the  fracture  surfaces  and  its 
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dependence  on  the  stress  level  was  evaluated  by  performing  fatigue  tests 
under  different  loads. 


v 


Table  /;  Loading  Conditions  for  Crack  Initiation  Times 
and  Damage  Evolution  Studies. 


*7  max 

(MPa) 

®  min  (MPa) 

Frequency 

(sec-1) 

Cycles  to 
Initiation  (xIO3) 

Cl: 

8.00 

2.75 

0.33 

300±1 0 

C2: 

10.50 

2.75 

0.33 

30±2.5 

C3: 

14.30 

2.75 

0.33 

4±0.5 

12.00 

3.00 

0.33 

20±2.0 

13.00 

3.00 

0.33 

1 0±1 .2 

3.  Results  and  Discussion 

Upon  application  of  the  load,  a  core  zone  of  crazes  appeared  at  the  notch  tip. 
While  this  zone  increased  in  size,  few  arched  crazes  nucleated  and  propagated 
around  the  core  zone.  As  time  progressed,  the  crazes  grew  longer  as  well  as 
new  crazes  appeared  in  and  around  the  vicinity  of  the  core.  A  series  of  optical 
micrographs  taken  during  a  typical  experiment  is  shown  in  Figure  2.  Crack 
initiation  was  defined  at  the  instance  when  a  crack  appeared  at  the  notch  tip. 

One  way  to  investigate  damage  evolution  is  to  compare  the  contours  of 
equal  damage  density  at  consecutive  configurations.  Such  pointwise 
comparison,  however,  may  not  be  appropriate  because  of  local  fluctuations  in 
damage  density.  Indeed,  the  contours  of  equal  damage  level  at  each 
configuration  obtained  from  different  specimens  exhibit  noticeable  fluctuation 
on  the  scale  of  several  micrometers  (Figures  3  to  7). 
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Cycles  38,760  Cycles:  39,430 


Figure  2  Series  of  optical  micrographs  showing  damage  growth  before 
crack  initiation. 
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Therefore  integrai  parameters,  such  as  the  moments  of  damage  distribution, 
were  used  here  to  evaluate  the  properties  of  damage  evolution  before  crack 
initiation. 

Characterization  of  any  distribution  with  the  use  of  moments  should 
involve  comparison  of  a  sufficiently  large  number  of  moments.  However, 
relatively  high  order  moments  may  not  be  appropriate  to  compare  because  the 
experimental  error  built  into  the  measurement  makes  them  relatively  inaccurate. 
In  the  present  analysis,  the  distributions  of  damage  were  characterized  by  the 
total  number  of  crazes  within  the  zone  lo,  the  second  kx.  l2y  and  fourth  Ux,  Uy, 
and  l2x2y  central  moments.  These  quantities  were  evaluated  from  the  following 
expressions, 

Ni  N2  Ni  N2 

l0=IIPij.  Imxny  =  £  £  [Xjj-Xc]m[yjj-yc]nPijAyjAxj 
i-1  j-1  i-1  j-1 

Here  i  =1  ....N-j  and  j  =1  ,...N2  refer  to  the  size  of  the  mesh  (see  Figure  1)  and  m, 
n  stand  for  the  order  of  the  moments  and  take  values  of,  0,  2,  or  4.  Xc  and  yc 
are  the  coordinates  of  the  damage  zone  center  with  reference  to  the  notch  tip 
and  given  by, 

Ni  N2  Ni  N2 

Xc  =£  X  XijPij  AVjAXi / X  £  Pij  AyjAXi ,  yc=0 
i-i  j-i  i-i  j-i 

yc  is  zero  due  to  symmetry  of  the  damage  zone  around  the  y  axis  (Figure  3). 

The  contours  of  equal  damage  density  p,  pertaining  to  two  specimens 
fatigued  under  conditions  C2  (Table  I)  and  up  to  the  instance  of  crack  initiation 
are  displayed  in  Figure  3.  Note  that,  the  contours  of  damage  density  have  been 
symmetrized  with  respect  to  the  bisector  of  the  V-notch  because  small  variations 
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in  density  were  observed  at  some  points  located  at  equal  distances  from  the  X 
axis.  These  differences  were  attributed  to  local  fluctuation  in  damage  density. 


y 
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Figure  3  Contours  of  equal  damage  density  p  [mm2/mm3],  for  two 

fcondjtfonsfC2)Ued  under  the  same  conditions  and  equal  times 
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Thus,  the  density  at  these  locations  was  substituted  with  the  arithmetic  mean  of 
the  respective  density  values.  For  the  sake  of  visual  clarity,  the  level  of 
experimental  error  is  not  shown  in  Figures  3  to  6.  The  error  in  the 
measurements  of  craze  density  far  from  the  notch  tip  was  about  10%.  For  the 
measurements  near  the  notch  tip,  the  error  was  about  15  -  20%.  The  larger 
value  of  error  was  because  the  high  craze  density  in  the  close  vicinity  of  the 
notch  tip  limits  optical  microscopy  which  makes  accurate  measurements 
difficult.  From  these  data  moments  of  the  distribution  with  respect  to  the  center 
of  the  zone  as  well  as  the  notch  tip  were  evaluated  (Table  II).  A  comparison  of 
the  moments  shown  in  Table  II  shows  that  the  largest  difference  is  -10%  and 
occurred  between  the  l4X  moments  with  reference  to  the  notch  tip.  The  rest  of 
the  differences  were  always  less  than  10%.  These  data  and  the  reproducibility 
of  damage  distribution  within  a  process  zone  in  the  same  material  [17]  were 
taken  as  sufficient  indications  that  damage  distributions  during  crack  initiation 
were  reproducible. 

Table  II:  Inertia  Moments  for  two  Configurations  of  Damage  before  Crack 
Initiation. 


Configuration 

1 

2 

1 

2 

Reference 

Notch  Tip 

Notch  Tip 

■nmsi 

55.0 

0.0 

52.0 

0.0 

l0  (number) 

732 

680 

lac  (mm2) 

2.8 

2.5 

0.65 

0.63 

Mmm2) 

0.46 

0.42 

0.46 

0.42 

Ux  (mm4) 

0.20x1 0'1 

0.1 8x1 01 

0.11x10*2 

0.11x10-2 

Ur  (mm4) 

0.57x10-3 

0.52x10-3 

0.57x10-3 

0.52x10-3 

'aav  (nun4) 

0.21x10-2 

0.18x10-2 

0.40x10-3 

0.39x10-3 
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Measurements  of  damage  density  within  the  zones  yielded  the  contours  of 
equal  density,  p,  that  are  shown  in  Figures  4,  5,  and  6,  for  conditions  Cl,  C2 
and  C3,  respectively.  From  the  experimental  data  presented  in  these  Figures 
the  center  of  the  damage  zone,  Xc,  can  be  calculated  at  each  configuration. 
The  data  in  Figure  7  indicate  that  the  center  of  gravity  of  the  damage  zone 
increased  monotonically  with  the  number  of  cycles.  To  investigate  the  type  of 
transformation  the  damage  density  is  undergoing  during  its  evolution,  the 
following  characteristic  scales  along  the  X  and  Y  axes  are  defined, 

=  ]2£  =  o§ 

lo  lo 

Since  1^.  l2y  are  the  second  moments  of  the  damage  zone  and  l0  is  the  zeroth 
moment  (which  expresses  the  total  amount  of  crazes  within  the  damage  zone) 
ox  and  ay  can  be  looked  upon  as  the  measurements  of  damage  spread  in  the  X 
and  Y  axes,  respectively.  Accordingly,  for  the  purpose  of  comparing  damage 
distributions  within  the  damage  zone  the  following  ratios  were  defined  between 
ith  and  jth  configurations  (i=  1 , 2, 3  and  j=i+1 ), 


The  type  of  transformation  of  damage  distribution  between  consecutive 
configurations  was  examined  by  comparing  the  following  quantities: 
(Xjj,  Aji)),(|Xjj,  Mjj)),  and  (A|JJ.ji,  Hji)).  The  above  defined  ratios  for  conditions 
Cl,  C2,  and  C3  are  shown  in  Table  III.  The  data  in  Table  III  indicate  a  relatively 
large  difference  in  the  ratios  between  the  first  and  second  configurations  for 
conditions  Cl . 
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Figure  7  Evolution  of  the  gravity  center  of  the  damage  zone. 
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This  may  be  because  the  system  needs  some  time  for  damage  to  be  developed 
so  that  statistical  measurements  are  meaningful.  However,  differences  between 
the  rest  of  the  ratios  were  small. 

Table  III:  Moment  Ratio  of  Damage  Distribution  at  Consecutive 
Configurations  before  Crack  Initiation. 


Condition  Cl 

EH99I 

E9BS9 

HU 

KHE9 

1HBI 

■199 

i 

msm 

B  9 1 

EBB8H9 

Condition  C2 

1S8R1IM 

BflSi 

BSII 

1 

■9911 

Condition  C3 

HI 

msm 

■RRE9I 

BBS 

H^i  *  2.10 

inn 

m 

HSHI 

The  constancy  of  the  ratios  Ap,  |Xji  between  configurations  shown  in  Table  III 
implied  that  the  evolution  of  damage  can  be  approximated  by  a  linear 

transformation  of  the  space  variables.  To  obtain  the  kinematic  parameters,  the 
velocity  Vm(x)  of  a  point  within  the  zone  is  expressed  by  the  first  two  terms  of 

Taylor  series  around  the  center  of  the  damage  zone, 

Vm(x)  ~  Vm(0)  +  Vm  n(x)xn  (1 ) 

Here,  the  first  term  represents  the  rigid  translation  of  the  damage  zone  and  the 
second  term  represents  the  rotation  (anti-symmetric  part  of  Vm  n)  and  the 
deformation  (symmetric  part  of  Vm  n)  which  in  turn  can  be  decomposed  into 
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isotropic  expansion  and  homogeneous  distortion  of  the  damage  zone,  in  this 
section  the  velocity  Vm(x)  is  expressed  in  terms  of  the  transformation  coefficients 
Xp  and  tip.  Subsequently  the  rates  of  expansion  and  distortion  are  evaluated. 

The  transformation  matrix  of  points  with  equal  damage  density  between 
configuration  i  and  j  is, 


D 


X$  0 

.0  m. 


(2) 


which  can  be  rewritten  as, 


D 


1  +AXp 
0 


0 

1+Apjj 


=  l+AD 


(3) 


where  I  is  the  identity  matrix  and, 


AD  = 


AA|  0 

0  Apf 


(4) 


Thus,  the  displacement  field  within  the  damage  zone  can  be  expressed  as. 


X+AX 

Ly+AyJ 


1  +AA,p 


0 

1+Apjj 


x 

y 


AXc 

0 


(5) 


where,  Axc,  represents  the  growth  of  the  center  along  the  x  axis.  Equation  (5) 
can  be  rewritten  as, 


Ax 

i 

> 

2? 

o 

J 

X 

4- 

’  AXC  ' 

.  Ay  . 

.  0  A*ijj  J 

y . 

V 

0 

Dividing  both  sides  by  the  increment  of  the  cycle  number  AN,  and  setting 

an ->Vl'^  ^V2'  Sr  S  * is  obtained’ 
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Vi 

V2 


Xjj  0 

0  lip 


(7) 


X  1 

y 


*c 

o 


A  comparison  of  Equations  (1)  and  (7)  results  in, 

Vm(x)  = 


[Vil 
1  1 

.  Vm(0)  = 

*c 

.  Vm,n(x)  = 

Xji  0 

lV2J 

L  0  J 

.0  Mji. 

(8) 


Thus,  the  rates  of  expansion  e,  and  distortion  d,  of  the  damage  zone  can  be 
expressed  as, 


e  -  and  d  - 


(9) 


The  average  values  of  e  and  d  as  well  as  the  growth  of  the  gravity  center  of  the 

damage  zone  for  three  loading  conditions  are  displayed  in  Table  IV. 

Table  IV:  Growth  Ratos  of  the  Gravity  Center,  Isotropic  Expansion  and 
Distortion  for  three  Loading  Conditions . 


Configuration 

AVAN  (x104) 
(jim/cyde) 

Ae/AN  (x106). 

1 /cycle 

Ad/AN  (xIO6), 

1 /cycle 

Conditions  Cl : 

1-2 

1.3 

13.0 

4.4 

2-3 

1.5 

3.0 

1.0 

3-4 

1.6 

2.0 

0.7 

Conditions  C2: 

1-2 

4.5 

20.0 

4.5 

2-3 

11.8 

15.0 

3.6 

Conditions  C3: 

1-2 

141.5 

27.0 

48.0 

2-3 

4.5 

1.75 

3.0 

The  data  in  Table  IV  show  that  while  the  growth  behavior  of  AXC/AN  depended 
upon  the  loading  history,  the  average  rates  of  the  damage  zone  deformation, 
Ae/AN  and  Ad/AN,  decreased  until  crack  initiation. 

It  is  worth  pointing  out  that  damage  evolution  within  a  process  zone  on  the 
same  material  can  be  described  by  a  linear  transformation  [17].  These  findings 
and  the  results  of  the  present  studies  indicate  that  a  similar  transformation 
relates  damage  evolution  before  crack  initiation  and  during  slow  crack  growth. 

The  evolution  of  an  average  damage  density  <p>  =  I/A,  where  I  is  the 
total  number  of  crazes  within  the  zone  and  A  is  the  area  of  the  zone  with 
normalized  cycle  number  is  shown  in  Figure  8.  (Note  that  because  of  the  large 
difference  in  the  cycles  to  crack  initiation  between  the  different  experimental 
conditions,  Cl :  300,000  cycles,  C2:  30,000,  C3: 4,000,  the  number  of  cycles  has 
been  normalized  with  the  respective  cycle  number  at  initiation  in  order  to  clearly 
see  the  trend  of  <p>).  It  is  interesting  to  note  that  all  three  sets  of  data  can  be 
approximated  by  an  equation  of  the  form  (lines  in  Figure  8), 

<p>  =  l-exp(jcN)  (10) 

where  N  is  the  number  of  cycles  and  K  is  a  load  dependent  parameter  and 
equal  to  K  =  2.5x1 0*5  (Cl),  8x1  O'5  (C2),  and  1.5x1  O'3  cycle'1  (C3).  Expression 
(1 0)  indicates  that  the  dependency  of  <p>  on  the  cycle  number  is  of  a  damping 
character.  Similar  trends  of  damage  accumulation  have  been  reported  on 
smooth  specimens  of  different  polymers  under  creep  loads  [20]. 

The  data  in  Figure  8  show  that  while  the  rates  of  damage  accumulation 
were  different,  craze  densities  <p>,  at  crack  initiation  were  40,000,  44,200  and 
47,000  mm'2,  respectively.  Note  that  the  error  between  the  extreme  values  of 
these  densities  was  16%  while  the  largest  difference  from  the  mean  was  about 
10%. 
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Thus  one  may  be  led  to  consider  the  value  of  <p>  as  being  independent  of  the 
loading  conditions  at  crack  initiation.  This  would  be  realistic  if  the  density  within 
the  zone  was  constant  throughout  the  zone.  However,  the  results  of  these 
studies  clearly  indicate  that  damage  was  not  homogeneously  distributed 
(Figures  3  to  6).  Thus,  an  average  quantity  like  <p>  may  not  be  indicative  of  the 
crack  initiation  event.  On  the  other  hand,  crack  initiation  occurred  within  a  core 
zone  ahead  of  the  notch  tip.  Morphologies  of  this  highly  localized  damage  zone 
at  initiation  for  conditions  Cl  and  C3  are  shown  in  Figures  9  and  10, 
respectively.  Within  the  resolution  of  the  experimental  measurements,  damage 
density  within  the  core  was  about  1 ,300  mm2/mm3  and  approximately  the  same 
for  all  loading  conditions  (Figures  4  to  6).  Similar  value  of  damage  density  has 
been  obtained  in  the  close  proximity  of  the  crack  tip  grown  under  low  frequency 
fatigue  loads  [17].  The  constancy  suggests  that  a  certain  level  of  damage 
density  is  required  for  crack  initiation  and  subsequent  crack  growth  and  that  this 
density  may  be  a  material  parameter.  Additional  experimental  as  well  as 
theoretical  research  is  needed  to  firmly  establish  this  important  observation. 

The  optical  micrographs  in  Figures  9  and  10  bring  out  an  interesting 
observation.  While  the  density  at  the  notch  tip  was  practically  the  same,  the 
cycle  numbers  to  crack  initiation  were  drastically  different,  i.  e.,  300,000  and 
4,000  cycles,  respectively.  Moreover,  damage  dissemination  around  the  core 
zone  was  quite  different  (Figures  9  and  10).  Thus,  one  would  expect  that  in  the 
case  where  large  crazes  surrounded  the  core  zone  would  take  longer  to  initiate 
a  crack  (Figure  10)  due  to  a  greater  extent  of  shielding  of  the  peripheral  crazes 
on  the  core.  However,  this  was  not  the  case.  Therefore,  it  is  assumed  that  time 
dependent  processes  within  the  damage  zone  play  a  more  important  role  than 
the  interaction  between  the  crazes. 
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Figure  9  Optical  micrographs  showing  damage  distribution  at  crack 
initiation  (Conditions  Cl). 
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Figure  1 0  Optical  micrographs  showing  damage  distribution  at  crack 
initiation  (Conditions  C3). 
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It  was  stated  earlier  that  a  small  crack  was  observed  at  the  notch  tip.  This 
event  was  associated  with  crack  initiation.  An  accurate  initial  crack  size, 
however,  was  difficult  to  detect  with  the  use  of  standard  optical  observations 
during  the  fatigue  tests.  To  evaluate  an  initial  crack  size,  the  specimens  were 
pulled  to  fracture  shortly  after  crack  initiation  and  the  fracture  surfaces  were 
observed  under  an  optical  microscope.  A  typical  micrograph  of  the  fracture 
surface  morphology  near  the  notch  tip  is  displayed  in  Figure  1 1 . 


Notch  Tip 


% 

* 

Note  that  two  distinct  morphologies  can  be  seen  in  the  fracture  surface:  (i)  a 
relatively  rough  surface  from  the  notch  tip  which  extents  up  to  an  arched  front, 
and  (ii)  the  surface  with  tearing  like  features.  It  is  assumed  that  the  first 
morphology  is  associated  with  crack  initiation.  The  second  morphology 
corresponds  to  the  crack  growth  phase  resulting  from  the  simple  pulling  of  the 
specimen  after  initiation.  The  morphology  shown  in  Figure  1 1  suggested  that 
the  crack  front  was  not  straight.  Instead,  crack  grows  more  in  the  middle  than  at 
the  edges  of  the  specimen  surface.  The  crack  length  at  initiation  { ,  taken  as  the 
average  of  five  measurements  from  the  notch  tip  to  the  curved  front  is  shown  in 
Figure  12  as  a  function  of  the  stress  amplitude.  These  data  indicate  that  the 
dependence  of  the  crack  length  on  the  applied  stress  is  of  the  form, 
Ln 4  =  Bq  +  Bia  where  Bo  and  Bi  are  constants. 


Figure  1 2  Crack  length  at  crack  initiation  for  various  levels  of  stress 
amplitude. 
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Crack  initiation  times  Tj,  plotted  against  the  stress  amplitude  are  shown  in 
Figure  13  (note  that  Nj  =  vt;  where  Nj  is  the  number  of  cycles  for  crack  initiation 
and  v  is  the  test  frequency).  These  data  show  that  Tj  can  be  approximated  by  an 
expression  of  the  form,  Lrttj  =  Do  -  Dig  where  Do,  and  Di  are  parameters. 
Similar  relationships  have  been  obtained  for  different  materials  including 
polymers  and  metals  under  constant  loads  and  different  temperatures  [21 ,22] 
and  have  been  explained  in  terms  of  stress  -  temperature  activated  processes 
at  the  molecular  level.  It  may  be  expected  that  such  a  relation  would  not  apply  in 
fatigue  loads.  However,  the  frequency  of  the  fatigue  experiments  was  very  low. 
Thus,  it  may  be  that  the  material's  response  is  primarily  due  to  the  time  under 
load  and  less  due  to  cycle  reversal.  Similar  results  were  obtained  under  cyclic 
loads  in  PANB  rocket  propellant  material  [23]  and  other  polymeric  solids  [24]. 


Figure  13  The  time  to  crack  initiation  plotted  against  stress  amplitude. 
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4.  Concluding  Remarks 

This  investigation  was  intended  to  study  the  phenomenon  of  crack  initiation 
under  low  cycle  fatigue  conditions  in  an  amorphous  polymer.  The  results  of  the 
work  outlined  in  this  paper  have  shown  that: 

1 .  A  damage  zone  accompanied  the  process  of  crack  initiation.  Within  the 
zone  two  patterns  of  damage  were  distinguished:  a  core  of  highly  dense 
damage  and  a  peripheral  less  dense  crazes  around  the  core  zone  Crack 
initiation  occurred  within  the  core  zone. 

2.  Damage  evolution  before  crack  initiation  evolved  in  a  self  -  similar  manner. 
That  is,  a  linear  transformation  related  the  points  of  equal  damage  at 
consecutive  configurations  of  damage. 

3.  While  the  growth  behavior  of  the  damage  zone  center  depended  upon  the 
loading  history,  the  average  rates  of  the  damage  zone  expansion  and  distortion 
decreased  until  crack  initiation.  Moreover,  the  dependency  of  an  average 
damage  density  with  the  cycle  number  exhibited  a  damping  behavior. 

4.  Damage  density  within  the  core  of  damage  was  found  independent  of  the 
loading  conditions,  however,  the  pattern  of  the  peripheral  crazes  was 
dependent  upon  the  loading  conditions.  The  large  difference  in  crack  initiation 
times  was  attributed  to  the  role  of  time  dependent  processes  within  the  damage 
zone  and  the  interaction  between  the  crazes. 

5.  The  crack  length  at  initiation  was  found  to  increase  exponentially  with  the 
applied  stress. 

6.  The  time  to  crack  initiation  was  related  to  the  stress  level  with  an 
exponentially  decaying  relationship.  The  results  agree  with  fracture  models 
based  on  reaction  theories. 

The  results  of  the  present  studies  and  those  reported  in  [17,  18]  indicate 
certain  similarities  between  the  characteristics  of  damage  before  crack  initiation 
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and  during  slow  crack  growth.  These  are,  (a)  damage  evolution  can  be 
described  by  a  linear  transformation,  (b)  the  densities  of  damage  within  a  core 
zone  at  crack  initiation  and  during  slow  crack  growth  are  of  the  same  level. 
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